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INTRODUCTION 

Work is currently in progress to obtain temperature distributions of dust in the most infrared-luminous 
galaxies. The results presented herein are of a preliminary nature, representing a zeroth-order approximation 
(sec WORK YET TO BE DONE, below). The objects which have been analyzed so far are all galaxies 
from the IRAS Bright Galaxy Sample with infrared luminosities Ljr > 10 u L e (see Carico el ai, 1988). 
They are: Arp 220, Mrk 231, Mrk 273, NGC 1614, NGC 3690, NGC 6285/6, and Zw 049.057. The 
analysis utilized 3.7 fim data from the Palomar 5 m Hale telescope, IRAS data at 12, 25, 60, and 100 /im, 
and 1 mm continuum data from the CalTech Submillimeter Observatory on Mauna Kea. 


METHOD 

A cloud of N spherical dust grains at a distance d from the observer produces an observed flux density 
at a wavelength A given by 

/„(A) = J nB v (X,T)(^j Qat.W — 0) 


where a is the radius of the dust grains, Q a 6»(A) is the absorption efficiency of the dust grain material, 
and B V {\,T) is the Planck function. If a mean mass density is adopted for the grain material (a value of 
3 g cm -3 was used to obtain the accompanying results), N can be rewritten in terms of the total mass of 
dust, M, and eqn. (1) can be put in the form 
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where C is a constant for a given emission source, depending only on the mean mass density of the grain 
material and the distance to the source. For A > a, the quantity Q a bs(X)/a is independent of a (see, e.g., 
Hildebrand, 1983); thus, for infrared emission from spherical interstellar dust grains, the observed flux 
density does not depend on the size of the grains. 

In numerical form, eqn. (2) becomes 


/i»(A<) = C 


QabAXi) 1 

a A 3 


£ 


, \ /kTj 



dTj . 


Using logarithmic intervals for A and T, so that A, = A 0 x', Tj = T 0 y 1 , and dTj -Tj — 7)_i - Tj( 1 — ^), 
where x, y, X 0 , and T 0 are constants, and defining 
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one obtains 


( 3 ) 
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which must be solved for the pj. 

For the analysis presented here, eqn. (3) was solved using a modified version of the least-squares 
computer program described by Pajot el al. (1986) (which was generously provided by J. L. Puget when 
this author’s own program failed to cooperate for as yet unknown reasons - my sincerest thanks!). A 
wavelength range from 3.7 /im to 1000 fim was used for all objects. The results are shown in the 
accompanying figures. 


RESULTS 

The results for three representative galaxies are shown in Figures 1 and 2. In each plot, three curves 
have been drawn, representing results for three different emissivity laws: 

Qabs OC A" 1 : solid line 
Qabs <x A -1 5 : dotted line 
Qabs oc A“ 2 : dot — dashed line 


where Q abs is the absorption efficiency, normalized to Q afc ,(100 fim)/a = 133.3 cm- 1 , and a is the 
effective grain size. 

Figure 1 shows the mass of dust (in units of Mq = 2 x 10 33 gm) per logarithmic temperature interval. 
The total mass is thus the area under the curve shown. 

In Figure 2 is shown the fraction of the total luminosity emitted by dust at a temperature T (in the top 
half of each plot), and the fraction of the luminosity which is being emitted by dust which is at temperatures 
> T. These plots tend to bring out the differences between the galaxies much more clearly than the mass 
distribution plots of Figure 1, due to the extreme temperature dependence of the luminosity (L oc T 4 ^, for 
Qabs oc \~P). In particular, one can compare the very sharply peaked distribution for Arp 220, centered at 
50 K, and falling to 5% of the luminosity at temperatures £ 25 K and £ 100 K, with the broad distribution 
for Mrk 231, which peaks at roughly 150 K and falls to 5% of the luminosity for temperatures ^ 40 K 
and £ 350 K. This difference can also be seen in the bottom plots where, for Arp 220, £ 10% of the 
luminosity is being emitted by dust with temperatures > 100 K, whereas for Mrk 231, the contribution to 
the luminosity from T > 100 K dust is ~ 60%. 


Figure 2 also indicates how uncertain estimates of the mass of cold dust can be. For the case of 
Mrk 231, Figure 2 suggests that more than 90% of the luminosity is being emitted by dust at temperatures 
£ 50 K. However, from Figure 1 it is seen that most of the estimated dust mass for Mrk 231 is at 
temperatures ^ 50 K. Hence, significant changes in the total mass of dust would not necessarily be 
reflected in the energy distribution. 


WORK YET TO BE DONE 

As mentioned previously, the work presented here is preliminary, intended to obtain a qualitative feel 
for the analysis and the results that can be obtained. The main simplifications which have gone into this 
analysis, and which will be addressed in subsequent work, are as follows: 

1. The dust was assumed to be everywhere optically thin to infrared radiation. This is almost certainly 
not the case, particularly for the compact central regions of such galaxies as Arp 220 and Mrk 23 1 . 

2. The 1 — 3 /im data has not yet been utilized, due to the complexities of accurately accounting for 
contamination from stellar emission. These wavelengths will clearly be important in understanding 
the hottest dust, particularly for temperatures in excess of 1000 K. 

3. The contribution to the emission from very small grains, presumably PAHs (polycyclic aromatic 
hydrocarbons), has not been addressed, and may be significant in some sources. 
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Figure 1: The mass of interstellar dust as a function of temperature for three infrared-luminous galaxies. The 
different curves on each plot are discussed in the text 
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luminosity emitted by dust at all temperatures > T. 
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